Abstract. The BL Lac Object 1ES 1426+428, at a red- 
shift of z=0.129, has been monitored by the Cat tele- 
scope from February 1998 to June 2000. The accumu- 
lation of 26 hours of observations shows a 7-ray signal 
of 321 events above 250 GeV at 5.2 standard deviations, 
determined using data analysis cuts adapted to a weak, 
steep-spectrum source. The source emission has an aver- 
age flux of $diff (400 GeV) = 6.73 ± 1.27 stat ± 1.45 syst x 
10 -11 cm -2 s _1 TeV _1 , and a very steep spectrum, with a 
differential spectral index of j—~ 3.60±0.57 which can be 
£Sj ■ refined to 7=— 3.66±0.41 using a higher flux data subset. 
<0 ' If, as expected from its broad-band properties, the Very 
, High Energy emission is hard at the source, these obser- 
vations support a strong absorption effect of 7-rays by the 
Intergalactic Infrared field. 
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1. Introduction 

BL Lac objects, together with Flat Spectrum Radio 
Quasars (FSRQs), constitute the extreme class of Active 
Galactic Nuclei known as blazars. They exhibit a highly 
variable non-thermal spectral energy distribution (SED) 
which is interpreted as the Doppler boosted emission of a 
jet pointing towards the Earth. 

While at least 66 blazars are repo rted in the third 
Egret catalogue ( Hartman et al. 1999 ), in the TeV range 
only two nearby BL Lac objects, Mkn 421 (z= 0.031) and 
Mkn 501 (z=0.034), have been firmly established so far. 
This can be explained partially by the small field of view of 
ground-based telescopes (a few degrees) , although, as com- 
pared to satellite-borne instruments, they benefit from a 
much larger sensitivity. The attenuation of 7-rays through 
pair-production with Intergalactic Infrared field (IIR) is 
another limiting factor for viewing distant sources in the 
Very High Energy (VHE) domain. However, the detection 
of absorption features in blazar spectra can provide inter- 
esting constraints on the poorly-measured 0.5 — 20 /im IR 
band. 

In this context, a survey of nearby BL Lacs 
has been carried out with the Cat telescope since 
1997. Among them, 1ES 1426+428 (a J20 oo=14 h 28 m 32.7 s , 
<5./2ooo=+42°40'21") at a red-shift of z = 0.129, occupies a 
peculiar position: the characteristics of its radio-to- X-ray 
SED are very similar to those of the two TeV BL Lacs, 
especially its hard X-ray emission which flags the presence 
of ultra-relativistic electrons. On the other hand, the rela- 
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tively high red-shift of 1ES 1426+428 can imply significant 
absorption by the IIR photons. 

VHE 7-ray detections of 1ES 1426+428 have recently 
been reported by Whipple and Hegra du ring 1999-2001 
( Horan et al. 2002 Aharonian et al. 2002] ). In this letter, 



we report on its detection based on observations made by 
the Cat telescope from February 1998 to June 2000. In 
Sect. 5 we will also give a first estimation of its spectrum. 

2. The CAT telescope and its standard data 
analysis 

The Cat imaging Cerenkov telescope, equipped with a 
3.1° high resolution camera (4.8° full field of view) and 
a 17.8 m 2 mirror, is located in the French Pyrenees at 



an altitude of 1650 m a.s.l. ( Barrau et al. 1998 ). The in- 
strument records Cerenkov light from the particles in 
cosmic-ray showers. The analysis is based on the compar- 
ison of individual events and theoretical average images 
of 7-rays as a function of impact parameter and energy 
(Lc Bohcc et al. 1998). To separate 7-ray initiated shower 
images from those due to hadronic showers the x 2 -like 
goodness-of-fit parameter, P x 2 , given by the above proce- 
dure is used, together with a cut on the pointing angle, 
a — which is the angle at the image barycentre between 
the actual source position and the reconstructed image 
axis of the 7-ray candidate — while requiring a minimum 
total charge, Q p . c in photo-electrons (p.e.). The standard 
cuts arc dedicated to strong and fairly hard sources ( "SH" 
hereafter). They are P x -i > 0.35, a < 6° and Q p . c > 30 p.e. 



(Piron et al. 2001), yielding a sensitivity of 4.5 standard 
deviations (<r) per -\/ Ion/ Hi for a Crab-like source at 
transit (~21° from Zenith) and a rate of 1.8 7 min -1 after 
cuts. 
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3. Weak point source data analysis 



Given the relatively high red-shift (£=0.129) of 
1ES 1426+428 as compared to the two confirmed 
BL Lacs, Mkn 421 and Mkn 501, one expects a rather 
weak flux and a soft spectrum due mainly to source 
distance and the IIR attenuation effect. The standard 
cuts have therefore been modified, trading-off 7-ray 
efficiency against background rejection. A study was 
made to get an indication of the IIR softening effect. 
Many estimates and models have been proposed for the 
IIR density (see Hauser & Dwek (2001) for a review). 



Conservative assumptions yield roughly a change of 
~ 1.0 — 1.5 of the spectral index for the red-shift con- 
sidered here. On the other hand, given its hard X-ray 
emission, a rather hard VHE spectrum is expected at 
the source, i.e. ~ —2. Consequently a weak simulated 
signal (1 7 min" 1 before cuts) with a spectral index ~ —3 
was tested against real background data taken from off 
observations of other sources (with roughly the same 
elevation as 1ES 1426+428). For such a steep spectrum, a 
tighter cut on P x 2 and a higher value of Q p . c were found 
to improve the background rejection, while a looser cut 
on a improves the efficiency because of the low-energy 
nature of the signal (for low energy 7-rays the estimation 
of the direction is less precise) . 
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Fig. 1. The accumulated ON-OFF pointing angle a plot 
for the three observing periods from 1998-2000. The OFF 
normalization has been done using the 20-100° control 
region. The increasing number of events at large a (>100°) 
corresponds to images cut by the camera border: the slight 
difference between ON and off samples in this region is a 
systematic edge effect which is not relevant to the signal 
or spectrum evaluations. 



Table 1. Cut values and results for MC generated signals 
for a source with a hard spectrum, after simulated IIR 
absorption. The first and second lines show results for the 
SH and WS cuts, respectively. e 7 is the 7-ray efficiency, 
Rh is the rejection factor for background cosmic rays, QF 
is the quality factor defined as QF = e 7 x \f~R~h- The 
final column gives the rate obtained on real Crab data in 
7 min -1 . Its variation from SH to WS is not the same as 
that of e 7 , due to the different spectral shape of the MC 
signal as compared to the Crab. 



Cuts 


P X 2 

min 


a 
max 


Qp.e 

min 


e 7 


R h 


QF 


7 min 1 
Crab 


SH 

WS 


0.35 
0.50 


6° 
8° 


30 
45 


32% 
25% 


180 
480 


4.3 
5.5 


1.8 
1.15 



The final values of the cuts ("WS" for Weak and Soft) 
were chosen by requiring a minimum efficiency of 25 % 
such as to keep a sufficient number of 7-ray events for 
spectral measurements and to avoid systematic effects due 
to a too small efficiency : P x 2>0.5, a<8° and <5 p . e >45 
p.e. The resulting quality factor, QF, and 7-ray rate on 
the Crab nebula, as compared to SH cuts, are given in 
Table [p. It has been verified that using the WS cuts for 
spectral analysis of Crab data yield, as expected, the same 
results as obtained by the standard cuts. 



4. Data sample and detected signal 

The observations reported in this paper were made with 
the Cat telescope from February 1998 to June 2000, in 
the ON-OFF mode, where control regions with the same 
trajectory on the sky are used as off-source data (i.e., with 
the same declination and shifted right ascension). 

A set of 26 hours ON and 12.5 hours OFF was selected 
according to criteria for good weather and smooth tele- 
scope operation (Piron 2000). Using the Cat data-base, 
the OFF sample was expanded by adding 59.9 hours of ob- 
servations of control regions (off) made during the same 
observing periods and with the same zenith angle range 
(i.e. 0-32°). The compatibility of the two OFF sample sets 
was checked over the tracking ratio distributions (i.e., the 
ratio of events within the a cut to the 20-100° control 



Table 2. Selected data and results for the three observing 
seasons from 1998 to 2000. The final column gives the 
significance of each detection according to the maximum 



likelihood definition by Li fc Ma (1983) 



Year 


ON (h) 


OFF (h) 


Excess 


7 min 1 


a 


1998 


8.9 


24.7 


91 


0.17 


3.1 


1999 


13.4 


30.5 


186 


0.25 


3.9 


2000 


3.7 


17.2 


44 


0.20 


1.8 


Total 


26.0 


72.4 


321 


0.21 


5.2 
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region) for all energy and zenith angle bins. This is the 
only relevant distribution to compare for signal and spec- 
trum calculations. Figure |l| shows the accumulated ON- 
OFF pointing angle a plot. An excess of 321 7-ray events, 
corresponding to a mean rate of 0.21 7 min -1 , is observed 
with a significance level of 5.2 a. The measured 7-ray rate 
is 18% of that of the Crab nebula within the same cuts 
(see Table |). There is a significant sensitivity improve- 
ment with the WS set as compared to SH cuts: the latter 
yield an excess of 367 events at 3.8 a. 

The detected signal in each observing season is sum- 
marized in Table ^. In 1998, an excess of 3.1c was seen in 
8.9 hours (February-May) with no flaring period evident. 

From February to June 1999, a 3.9a signal was de- 
tected in 13.4 hours. This period was consecutive to the 
campaign by BeppoSAX in February 1999 where the X- 
ray flux was very low as compared to archive data, but 
showed an extremely hard X-ray spectrum with a differen- 
tial index of —0.92 (Costamante et al. 2001). The source 



flux was the highest during two nights in March and April 
with some indication of rapid variability (see Fig. [2]): a 
rate of 0.8 Crab (3.4cr) followed by 0.3 Crab (1.5a) was 
detected in two runs (each of 30 min) at MJD 51261.5 
(23-24/03/99). During MJD 51278.5 (9-10/04/99) a rate 
of 0.6 Crab {2.7a) was followed immediately by two runs 
showing no signal. 
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Fig. 2. Event rate of 1ES 1426+428 during March and 
April 1999 in 7 min -1 . The source flux was the highest 
during MJD 51261.5 (0.8 Crab, 3.4a) and MJD 51278.5 
(0.6 Crab, 2.7a). The dashed line shows the average emis- 
sion level. 



Weather conditions in the Pyrenees limited the obser- 
vations in the ON - off mode to only 3.7 hours during 
spring 2000 (some other data are available in 2000, but 
with a different observation mode which is not analysed 



here) . The measured 7-ray rate is nevertheless compatible 
with that of the preceding periods. 

5. Spectral measurements 

Cat's spectrum analysis method has been applied 
to the data. The method (Djannati-Atai et al. 1999, 



Piron et al. 2001) relies on a maximum-likelihood esti- 



mation of the parameters taking into account the Pois- 
sonian probability distributions of the 7-ray and back- 
ground events, as well as the instrument response func- 
tions (effective area and resolution) derived by Monte- 
Carlo simulations and validated by comparison against 
nearly pure 7-ray samples ( Piron 1999[ ). Data are clas- 
sified within zenith angle and energy bins. The high en- 
ergy bins show no excess: when requiring the energy per 
event to be less than 1.3 TeV the significance level rises 
to 6.1a while almost no variation in the total number of 
excess events is seen. Consequently, a power-law is ad- 
justed in the limited range of 250 GeV to 1 TeV, yield- 
ing a differential flux O (400 GeV) = 6.73 ± 1.27 stat in 
10 -11 cm -2 s -1 TeV -1 units, and a very steep differential 
index 7 = —3.60+0.57. With no cut on the maximum en- 
ergy (i.e., from 250 GeV to 10 TeV) the fitted value of the 
index increases: 7 = —4.04+ 0.52. The corresponding in- 
tegral flux has been computed above the energy for which 
the errors on the two parameters 4>o and 7 de-correlate: 



5 int(>300 GoV) 



2.17+0.39 x 10 -11 cm- 2 s- 1 . 



To get a better determination of the spectral index, 
one can improve the signal-to-background ratio by select- 
ing the data, requiring a minimum 7-ray rate (High Flux 
subset). We make here the reasonable assumption that 
the background spectrum variations as a function of the 
statistical fluctuations of the tracking ratio are negligible. 
On the other hand, if the intensity-hardness correlations 



seen on other BL Lac sources (Djannati-Ata'i et al. 1999 



Aharonian et al. 2001) hold for 1ES 1426+428, the selec- 



tion could in principle bias the time-averaged spectrum. 
However, given the weakness of the signal and the magni- 
tude of a possible spectral variability (< 0.5) as compared 
to that of the statistical error, the bias should be negli- 
gible here. In consequence, a subset of 10.7 hours of ON 
data with a minimal 7-ray rate of 0.2 7 min -1 was selected 
(corresponding to an average rate of 0.5 7 min - ), and a 
differential spectral index 7 = 3.66 ± 0.41 was derived. 
Figure ^ shows the 68% confidence-level contour given by 
the likelihood estimator as well as the residuals to the fit 
for this subset. 

Systematic errors could arise from the uncertainty 
on the absolute energy scale (variations of the atmo- 
sphere transparency and light-collection efficiencies dur- 
ing the observation periods) and, to a lesser extent, 
from limited Monte-Carlo statistics in the determination 
of the response functions of the instrumen t. They hav e 
been estimated from detailed simulations ( Piron 2000 ): 
(A<V<M syst = ±20% and (A 7 ) s y st = ±0.06. Here, the 
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Fig. 3. 1ES 1426+428 time-averaged spectra between 
0.25 and ITeV from 1998 to 2000, for a power-law hypoth- 
esis using a High Flux subset with a mean 7-ray rate of 0.5 
7 min -1 . The areas show the 68% confidence-level contour 
given by the likelihood fit: O = 13.85 + 1. 77 stat ±2. 77 syst x 
(£/400GeV) 



-3.66±0.41" 



in 10~ n cm- 2 s _1 TeV- 1 units. 



The lower plot gives the ratio, in each bin of estimated 
energy (as opposed to true energy bins in the upper plot), 
of the predicted number of 7-ray events to the observed 
signal. The energy at which the errors on flux and index 
de-correlate is 430 GeV. 

systematic error in the index remains negligible compared 
to the statistical error, while for the flux these errors are 
comparable in magnitude. 



6. Discussion 

Observations of 1ES 1426+428 from 1998 to 2000 by Cat 
show a signal at 5.2er with an average event rate of order 
of ~ 0.2 Crab. This is the most distant detected BL Lac 
in the VHE range. 

For the three observing seasons studied here the source 
average emission seems stable, within the measurements' 
accuracy, while there is some evidence for time variability 
with transient emission rates > 0.5 Crab. 

Given its broad-band properties, 1ES 1426+428 was 
one of the most promising VHE candidates, despite its 
red-shift of > 0.1 ( Costamante et al. 2000p. Its detection 



supports the unifying scheme by Ghiscllini et al. (1998) 
where the lower luminosity blazars with high frequency 
peaked synchrotron emission (HBLs) are efficient accelera- 
tors to very high energies. It is remarkable that X-ray mea- 
surements flCostamantc et al. 2001 ), contemporaneous to 
Cat observations in February 1999 reported a peak emis- 
sion energy of ~ 100 keV, comparable to that of Mkn 501 
during its strong activity in 1997. 



If, as suggested by these data, the same acceler- 
ation/cooling mechanisms are at play within the two 
sources under similar conditions, one would expect a hard 
VHE differential spectrum with 7 ~ 2.0 at the source 



(see Djannati-Atai' et al. 1999). The very steep spectrum 
observed here, 7 = 3.66 ± 0.41, supports a strong ab- 
sorption of 7-rays by the diffuse Intergalactic Infrared 
field, even though the effect is compatible with several 
— among many — estimates and models of its density in 
the 0.5 — 20 fim band. Any definite conclusion on this sub- 
ject and any measurement of the IIR density will require 
the detection of more sources at different red-shifts. 
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